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We present a model describing the use of ultra-short strong pulses to populate the excited level of 
a two-level quantum system. In particular, we study an off-resonance excitation with a few cycles 
pulse which presents a smooth phase jump i.e. a change of the pulse's phase which is not step-like, but 
happens over a finite time interval. A numerical solution is given for the time-dependent probability 
amplitude of the excited level. The enhancement of the excited level's population is optimized with 
respect to the shape of the phase transient, and to other parameters of the excitation pulse. 
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I. INTRODUCTION 

Ultra-strong pulses with intensities of the order of 
10 15 W/cm 2 , and duration of the order of attoseconds, 
with just few optical cycles, are feasible with present day 
technology (see e.g. QrO). This technological develop- 
ment has been motivated by the large number of possible 
applications, several of which rely on coherent popula- 
tion transfer techniques. A partial list of such applica- 
tions is: stimulated Raman adiabatic passage (STI-RAP) 
J6H21/ adiabatic rapid passage (ARP) 1 10 ], Raman chirped 
adiabatic passage (RCAP) ITTI IT2l , temporal coherent 
control (TCC) JT3J HU, coherent population trapping 
lTT5l li"6l , optical control of chemical reactions [171 [18] , 
electromagnetically induced transparency (EIT) 1T51IT91 - 
22 1, efficient generation of XUV radiation 1 23 -26 ], coher- 
ent Raman umklappscattering|27J, breakdown of dipole 
blockade obtained driving atoms by phase-jump pulses 
1281 . Moreover, recently two schemes for efficient and 
fast coherent population transfer have been presented 
|29l , which use chirped and non-chirped few-cycles laser 
pulses. Another recent application [30] presents high- 
order harmonic generation obtained with laser pulses 
with a 7T-phase jump. Finally, the field of quantum in- 
formation processing benefits from these results, since 
many qubit realizations rely on precise quantum levels 
manipulation 131H34I . 

The presence of few optical cycles in the pulse gives a 
constant phase difference between the carrier wave and 
the pulse shaped envelope [35], in contrast with many 
cycle pulses [28 . 36 , 37 ]. Moreover, optimizing the pulses 
parameters is proven to enhance the excited state pop- 
ulation [38 1 or optimizing coherence in two-level sys- 
tems (TLSs) 1391 . In previous works we have already 
presented an analytical solution for the dynamics of a 
TLS excited with pulses of arbitrary shape and polariza- 
tion [40, 41]. Since in the model we present the change 
rate of levels' populations within a single optical cycle 



is not negligible, the rotating-wave approximation can't 
be used. In other words in the present model we can't 
neglect the contribution of the counter-rotating terms in 
the Hamiltonian ||40ll4H . 

In a previous work [42] we have presented a similar 
model, representing the interaction of a TLS with few- 
cycle pulses, where at time t = to the phase of the carrier 
wave jumps of an amount <ft, this jump being sharp and 
step-like. In that work, numerical analysis of the ana- 
lytic model has lead to an enhancement by a factor of 
10 6 - 10 8 in the population transfer, with the optimal 
phase jump of <fi — Tl and the optimal time coincident 
with the peak of the envelope. In the present work we 
improve that model, considering a smooth phase change, 
i.e. not step-like but happening over a finite interval of 
time. This new model more closely describes a realistic 
experimental scenario. 

The pulse is characterized by: Rabi frequency Qo, 
pulse width t, carrier frequency v, phase jump ampli- 
tude <p, phase jump time fo and phase jump duration 
At. Moreover, we consider two qualitative parameters: 
the phase jump shape, and the pulse envelope shape. 




Yt 

FIG. 1: The three functional shapes of the smooth phase 
jumps used: the red is a dropping hyperbolic tangent: <p(t) = 
(n/2)[l - tanh(5ai)], the blue is a rising hyperbolic tangent: 
<p(t) = (7i/2)[l - tanh(5af)], and the black one is a hyperbolic 
secant: (p(t) = {nil) sech(at)- m all the simulations the numer- 
ical normalized value is a = 0.265 



'Email: pkjha@physics.tamu.edu 



2 



EG 

o 

O 

W 



— j— i— i— i— i— j— i— i— i— i— j— i— i— i— i— r 




FIG. 2: Example of an excitation field with few oscillations and 
two different phases, thick solid (cosine) and dashed lines(sine), 
encompassed by a pulsed envelope (thick dashed line). Here 
the field is written in the form E(t) = £(f) cos(vf + cp ) where 
cp = for cosine and cp = -n/2 for sine pulse. 



We present an analytical solution for the time evolution 
of the excited state's population, together with a nu- 
merical simulation. In the numerical simulation we use 
3 functional shapes for the smooth phase jump: rising 
hyperbolic tangent, dropping hyperbolic tangent, and 
gaussian peek, (see figure 4), whereas for the envelope 
a gaussian peak has been used. Numerically optimizing 
the pulses parameters we have obtained enhancements 
for the population transfer of the order of 10 4 . 



A. NUMERICAL SIMULATION AND DISCUSSION 

Be \a) and \b) the states of a two-level atom (TLA), 
with energy difference Hco, and atomic dipole moment 
p. If we let this system interact with a classic field 
E(f) = £(r)cosvf, the equations of motion for the rela- 
tive wavefunctions are |43|: 



C a = i^l cos (vt)e^C b , 
n 

n 



(la) 
(lb) 



where A = ca - v is the detuning from resonance. Sim- 
ilarly to E3, defining /(£) = C a (t)/C b (t) and Q(f) = 
pS(t)/H, we have the following Riccati equation: 



/ + ;Q*(f)cos(vf)e" m 7 - iO.{t)cos(vt)e lan = 0. (2) 

The approximate solution for Eq. |2}, in terms of the tip 
angle 9 is given as in ] 40 1 
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where the tip angle 0(f) has been defined as 
0(f) = f D.(t')cos(yt')e lat ' dt' 

*J — DO 



(4) 



from which we have |C fl (f)l = |/(f)|/ Vl + 1/(0 1 2 - What is 
of interest is the asymptotic behavior of |C fl (oo)|. In 11421 
is shown good agreement between the analytical and a 
numerical simulation. To introduce the phase jump, we 
can write the Rabi frequency as 



Q(f) = Q (£)cosvf e^e'f 



ia,tM(t) 



(5) 



and then, using the same method as in 11421 , we can 
obtain an approximated analytic solution for the Riccati 
equation Q: 



fit) 



= i f df'<D(f')exp 2 f Q(t")dt" 



(6) 



The approximate analytical solution is in good agree- 
ment with the numerical simulation obtained by directly 
solving the coupled differential Eq.(l). From Fig. (3) we 
see that even for complex phase function the agreement 
is good. For the sake of completeness, we have added 
an appendix in which we show the strength of this ap- 
proach beyond standard TLS. Indeed the Riccati equa- 
tion approach gives a closed compact form from which 
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FIG. 3: Population left on the upper-level |fl)(b,d,f) as a func- 
tion of ratio of the carrier-frequency(v) of the excitation pulse 
and atomic transition frequency(£y c ), in the long time limit 
t » x, for corresponding phase jump function (a,c,e). Here 
the dashed line if the numerical simulation of Eq.(l) and the 
solid line is the approximate solution given by Eq.(6). For the 
excitation pulse we have used the form Qo(0 = Ae"" 2 ^ e l W> , 
and the phase functions have the following forms: (a) <p(t) = 
(n/2) sech[at], (c) <p(t) = n/2(sech[10af] + (1 + tanh[af]), (e) 
<p(t) = 7i/2(sech[af] + (1 - tanh[10af])). For numerical simula- 
tions we chose A = 0.035<i), a = 0.265y and y = 1.25a; where 
cj = (2n) 80 Ghz 
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FIG. 4: In this figure we present the results of the numerical analysis. Each of the three rows of plots refers to a different functional 
shapes of the smooth phase jump (phase change function). For each row we have a plot of the smooth phase jump, a plot of the 
excited state's population as function of time, and a plot of the excited state's population left after the pulse is gone as a function 
of the normalized excitation frequency. Similarly to figure|3] the functional shape of the excitation pulse is O (t) = Ae"" 2f V^ (f) . 
Moreover, for the plots of the excited state's population as function of time we have used the numerical value ofv/oj = 0.75. Phase 
change of the form (a) <p(f) = (n/2)[l + tanh(aii)], with three different values of a.\: red (steeper) a.\ = 5a; blue (in-between) ct\ = a; 
black (smoother) a.\ = 0.5a. (b) Corresponding behavior of the excited level population |C„(f)l- (c) Asymptotic value of the excited 
state population, as a function of the "resonance ratio" (excitation's frequency divided by transition's frequency) for this form of 
the phase change, (d) Phase change of the form cp(t) = (n/2)[l - tanh(«i f )], with (as in (a)) three different values of a\ : red (steeper) 
«i = 5a; blue (in-between) a\ = a; black (smoother) a\ = 0.5a. (e) Corresponding behaviour of the excited level population |C„(f)l- 
(f) Asymptotic excited population for this form of the phase change, (g) Phase change of the form: (p(t) = (7i/2)sech 2 (a 1 t), a\. 
red (larger) a\ = a; blue (in-between) a\ = 10a; black (narrower) a\ = 20a. (h) Corresponding behaviour of the excited level 
population |C a (f)|. (i) Asymptotic excited population for this form of the phase change. The value used for a is a = 0.265. For 
numerical simulations we chose A = 0.04375w, a = 0.265}' and y = 1.25cu where a> = (2n) 80 Ghz 



both the temporal and steady-state behavior of the two 
and three-level system can obtained. 

An interesting observation is that it is possible to 
rewrite the Rabi frequency in |5]l as 

Q(£) = Q (f) cos vte ,lc " t+ ^ t)] (7) 
and then define a>(t) = co + <p(t) / 1 and interpret this as a 



modulation of the atomic frequency, instead of a modula- 
tion of the excitation. Experimentally this can be realized 
in several ways, e.g. using modulated Zeeman or Stark 
effect. 

Now we move to discuss our numerical simulation 
of the dynamics of the two-level atom interacting with 
ultra-short, off-resonant and gradually changing phase 
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<p(t). We have performed numerical solution of the 
Riccati equation, using different types of phase change 
(smooth phase jump) functions. The result of this nu- 
merical analysis is shown in figure [4] The goal of this 
study is to find the best phase change which allows for 
the best coupling (most efficient energy exchange) of the 
excitation pulse with the excited state. 

In figure |4] we present the results of the numerical 
analysis. Each of the three rows of plots refers to a 
different functional shapes of the smooth phase jump 
(phase change function). For each row we have a plot 
of the smooth phase jump, a plot of the excited state's 
population as function of time, and a plot of the excited 
state's population left after the application of the pulse as 
a function of the normalized excitation frequency. Sim- 
ilarly to figure |3j the functional form of the excitation 
pulse is Clo(t) = Ae~ a2 ' 2 e'^ . Moreover, for the plots of 
of the excited state's population as function of time we 
have used the numerical value of v/co = 0.75. 

For this numerical simulation we have considered 
the following three phase functions (a)(b)(c): cp(t) = 
(tt/2)[1 + tanh(af)L (d)(e)(f): cj)(t) = (tt/2)[1 - tanh(at)] 
and (g)(h)(i) cp(t) = (n/2) sech 2 (af). We can see how the 
phase change duration At, i.e. the steepness of the (p(t) 
function, has not an unique effect on the excited popu- 
lation, and depends on the general shape of the phase 
change. In particular, it is worth noting that for ascend- 
ing and descending phase changes built on the tanh(£) 
function the effect of the steepness is opposite. We can 
observe a global behaviour which relates the characteriz- 
ing parameters of the phase change with the amplitude 
of the population of the excited state. Qualitatively, for 
the ascendent hyperbolic tangent we observe that by 
increasing the slope the population increases. On the 
other hand, for the descendent hyperbolic tangent the 
effect of this parameter is reversed: decreasing the slope 
of phase change leads to a decrease of the population. 
We remark that these behaviors are only global, and are 
reversed for some small ranges of frequencies. As an 
example, in plot 4.(f), for low ranges of laser frequen- 
cies, by decreasing the slope we increase the population, 
which is opposite of the behavior observed for higher 
frequencies. For the peaked shape (sech(£)), no general 
behaviors are observed. However, for the intermediate 
range of frequencies it can be observed a link between 
the increasing of the pulse width and the increase of the 
population. 

II. CONCLUSION 

To conclude, here we report our analytical and nu- 
merical results to show the effect of smooth phase jump 
on the dynamics (both transient and steady-state). We 
observed that the temporal profile of the phase jump 
function <p(t) has a profound effect on the excited state 
population \a). The two-level considered here can be 
the Zeeman-sublevels and the ultra-short (few to multi- 



cycle) pulse would be in the radio-frequency regime 
which has been reported in|44-46|. For an optimized 
phase function (within the set of parameters considered 
here), we were able to an observe an enhancement of 
10 4 in the population transfer. Such enhancement is 
seen in Fig. 4(f) (blue curve) which takes the value 
of |Cool = 0.002251 at v/co = 0.5 and phase function 
is (p(t) = (h/2)[1 - tanh(trf)] where a = 0.265. When 
ot\ = 5a = 1.325 the excited state population is enhanced 
by factor of ~ 2 X 10 2 . Similar enhancement is also ob- 
served for the phase function <p(t) = (n/2)[l - tanh(a£)] 
at frequency v/co ~ 0.9(see Fig. 4c). We not only can 
enhance excitation but for the same phase function and 
other choice of the parameter a, we can also suppress 
it. For example at near resonant excitation (see Fig. 
4(f)), the excited state population can be suppressed by 
~ 15-fold when a - 1.325 (red curve) in comparison 
to a.\ — 0.5a (black curve). Such control over excited 
state dynamics using smooth phase jump as an exter- 
nal parameter can be useful in microwave controlled 
Raman[47, 48], EIT with superstructures[49j to name a 
few while on the other hand with recent proposal on 
coherence-enhanced spaser|50 1 and propagating surface 
plasmon polaritons) 51 1 also viable areas to explore phase 
effects. The approximate analytical solution are in ex- 
cellent agreement for both delta function[42j or smooth 
phase jump considered here. We also extended this ap- 
proach beyond two-level atom to three-level in lambda 
configuration. 
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Appendix A: Analytical solution for three-level atom 

Motivation to add an appendix on the approximate 
analytical soltuion for three-level system in lambda con- 
figuration is to enlighten the strength of the method used 
to find the solution for two level atom with and without 
phase jumps. For the sake of simplicity we will con- 
sider constant phase <p. Let us consider a three-level 
atom(ThLA) in A configuration [see Fig. 5 inset]. The 
transition a <-> c is driven by the field Q2, while the field 
Qi couples the a o b transition. For the time scale con- 
sidered in this problem, we have neglected any decays 
(radiative and non-radiative). The equation of motion 
for the probability amplitudes for the states \a), \b) and 
\c) of the ThLA can be written as 

C a (t) = ihx{t)C b {t) + iCl 2 (t)C c (t) (Al) 
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FIG. 5: Numerical(red dotted line) and analytical (blue solid 
line) solutions of the amplitude of the state \a) after long time 
in function of v/co for the laser pulse envelopes Qi(f) = Q2W = 
Q sech(at). For numerical simulation we chose Q = .04a;, a = 

0.075&), (Jj„b = <jO ac = CO = 1. 



q,(f) = iCi{(t)c a (t) 



(A2) 



C c (t) = iti* 2 (t)C a (t) (A3) 
where Q;(f) is defined as the effective Rabi frequencies 

&j(t) = Q ; <f) cos(v ; i)e !£U ''; j = 1,2 (A4) 
To solve for C„(f) and C c (f) let us define 



c b (ty 



C b (t) 



(A5) 



In terms of /(f) and g(t) Eqs.(|AT][A2][A3} reduces to 

/'(f) + iCTJ 2 (t) = Cl 1 + iClQ 2 g(t) (A6) 



g(t) + iCl\f(t)g(t) = iti* 2 f(t) 



(A7) 



In order to solve these equations we extended the 
method developed|40 41 ]. By neglect ing the non-linear 
term f 2 (t) and the term oc g(t) in Eq.l A6 i we can solve 
for /i(f) as 



f 

%J -co 



(A8) 



Similarly by neglecting the term oc g(t) in Eq.| |A7| we can 
solve for g\(t) as 



gl (t) = - f QWmt'W 

%J —oa 



(A9) 



where the tip angle 6\{t) is defined as 

di{t')= f Qxit'W (A10) 

U -co 

Next let us write the non-linear term in Eq.(|A6) as 



f(t) = \f(t) - m)f + m)M) - fi(t) (ah) 



Then Eq.| |A6) can be written as 

f(t) + iCl\(t){[f(t) - MO? + 2/(f)/i(f) - flit)} 

= jQ a (f) + i€l 2 (t)g(t) 



(A12) 



Let us assume that g(t) » gi(t) and w e neg lect [f(t) - 
/i(t)] 2 |40| in this case we can write Eq.(A12 1 in term of 
the tip angles 6i(t) and 6 2 (t) 

f(t)+i&[®{2f(t)Mt)-Ji®} = iUt)+iUt)gi{t) (A13) 
where 



e2in= r 

•J — oc 



Cl 2 (t')dt', 



(A14) 



The analytical solution of the equation Eq.l A13 1 is then 



/(0 



~ a(t) f b{t')e a(t,) dt' 

Jf 



where 



a(x)=2i9 1 (t)f 1 (t) 



and 



b(x) = + ie 2 (t) gl (t) + i&[(t)ft(t) 



(A15) 



(A16) 



(A17) 



For g(t) the solution can be obtain from Eq.|A7lwhere 
we use /(f) « fi(t) : 



g(t) + ie\{t)h{t)g{t) = id 2 (t)Mt) 

which give us 



git) 



-c(t) 



f D(t')e c{t,) dt' 



where 



and 



c(x) = i&[(t)Mt) 



D(x) = W 2 (t)Mt) 



(A18) 



(A19) 



(A20) 



(A21) 



In Fig. 5 we have plotted the Numerical(red dotted 
line) and analytical (blue solid line) solutions of the 
amplitude of the state \a) after long time in function 
of vjco c for the laser pulse envelopes Qi(f) = 0. 2 (t) - 
□osech(af) with Qo = .04oj,a = 0.075co,a> a b - co ac - 
co = 1. We see that the approximate analytical so- 
lution matches well with the numerics under the pa- 
rameters considered here. Extension of this methodol- 
ogy to Schrodinger equation[52 , 53] and position depen- 
dent mass Schrodinger (PDMSE) equation can be found 
in ||54ll55l . 
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